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ABSTRACT 

We consider the cosmology of modified gravity models in which Newton's 
constant is distorted by a function of the inverse d'Alembertian acting on the 
Ricci scalar. We derive a technique for choosing the distortion function so as 
to fit an arbitrary expansion history. This technique is applied numerically 
to the case of ACDM cosmology, and the result agrees well with a simple 
hyperbolic tangent. 
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1 Introduction 



Evidences in favor of an accelerating cosmic expansion are now quite numer- 
ous, reaching from the first systematic observations of Type la supernovae [1] 
to the more recent WMAP survey of the Cosmic Microwave Background [2] . 
The standard interpretation of this acceleration, in the form of a cosmological 
constant [3], raises major well known questions. In particular it is not under- 
stood why the observed value of the vacuum energy density associated with 
the cosmological constant is so close to the non relativistic matter energy 
density (the so called coincidence problem) and so far from the "natural" 
order of magnitude expected from high energy physics [I]. 

The problem is that the geometry reconstructed from observation is not 
sourced by known matter according to the Einstein equation, 

(^-i^fl) obs ^(8.G N ^) t _. (1) 

Different ways addressing this problem can be classified as to whether it 
is the right hand (matter) side or the left hand (gravity) side of ([1]) which 
is modified. "Dark energy" models explain the data by introducing a new 
source of stress energy to the right hand side of (pQ). A cosmological constant 
is one example; another is a scalar quintessence field ip described by the 
Lagrangian [5], 

£ = - l -d^d v ^^Tg - V{<p)yf^ . (2) 

Local, metric-based modifications to the left hand side of (CQ) are restricted 
by stability |6J to take the form of replacing the Ricci scalar of the Hilbert 
Lagrangian by an arbitrary function of the Ricci scalar [7J [8] , 

~ TTr- R ^ = 9 ^ -TTr- F ^y/9 • (3) 

More general modifications of gravity must either abandon stability or local- 
ity [91 [Till EEH [12] , or they must involve some field other than the metric to 
carry part of the gravitational force [131 EE]- Other proposals involving ex- 
tra dimensions or massive gravitons have not yet reached a fully satisfactory 
state [T51IT5]. 

Central to the evaluation of any class of models is the reconstruction 
problem. This consists of identifying the extent to which parameters in the 
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model such as V(ip) and F(R) can be adjusted to support a geometry of 
the Friedmann-Lemaitre-Robertson- Walker (FLRW) form with flat spatial 
sections, 

ds 2 = -dt 2 + a 2 (t)dx-dx , (4) 

where the scale factor a(t) is a known but arbitrary function of time. Of 
course there is just one expansion history in nature and a putative model 
need only explain that, but if reconstruction can be solved for a general scale 
factor then it is certain a model within the given class can fit the actual 
expansion history. To the extent that the solution is constructive one also 
obtains important constraints on the model, of course limited by the precision 
with which the expansion history can be measured. 

It is straightforward to solve the reconstruction problem for scalar quint- 
essence models [HI 113 E] and a brief presentation of the solution will help 
focus ideas. For the FLRW geometry (J4]) the scalar must be independent of 
space, and only two of Einstein's equations are nontrivial, 



3H 2 
-2H - 3H 2 



8ttG n 
8nG N 



-V 2 + V((p\ 
\p 2 ~ Vltp) 



(5) 
(6) 



Here and henceforth, Gn is the Newton constant, a dot means a derivative 
with respect to the cosmic time t and we shall henceforth employ the usual 
definition of the Hubble parameter, H = a/a. By adding (jSJ) and (jSJ) one 
obtains the relation, 

- 2H = 8nG N cp 2 . (7) 

Hence one can reconstruct the scalar's evolution, and even invert it to express 
time in terms of the scalar, provided the Hubble parameter is monotonically 
decreasing, 



<p(t) 



tp ± / du\ 
Jo 



-2H(u) 



\ 87TG 



N 



t = t(<p) . 



(8) 



One then determines the potential by subtracting ([6]) from (jSJ) and evaluating 
the resulting (assumed known) function of time at t(<p), 



2H(t) + 6H 2 (t) 



16irG N V(<p) . 



(9) 
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Similar reconstruction procedures exist for F(R) models [HI IS] - The pur- 
pose of this paper is to solve the reconstruction problem for a recently pro- 
posed nonlocal cosmology model in which one multiplies the Hilbert La- 
grangian by an algebraic function of the inverse scalar d'Alembertian acting 
on the Ricci scalar [12J, 



The motivation for this class of models is to trigger late time acceleration 
by the transition from radiation domination, during which the Ricci scalar is 
nearly zero, to matter domination at about 10 5 years after the Big Bang. The 
subsequent time lag to the observed onset of acceleration, at about 10 9 years, 
would be provided by the effect of the transition being reflected through the 
nonlocal, inverse d'Alembertian. Although there is hope of eventually deriv- 
ing a model of the class fflQj) from quantum field theoretic loop corrections, 
the proposal is at this stage purely phenomenological. 

The free parameter f(X) in expression (fit)]) is known as the nonlocal 
distortion function. Absent a derivation from fundamental theory, it has the 
same status as the potential V(<p) in (j2J) and the function F(R) in (jHJ). What 
we will do in this paper is first to show how f(X) can be tuned to give 
an arbitrary a(t), then we will work out the specific form f(X) must take 
to reproduce the a(t) of ACDM, without actually employing a cosmological 
constant. This problem has already been studied in an excellent paper by 
Koivisto [IB] (see also [TH]) but for a local variant of the model, introduced 
by Nojiri and Odintsov [2U], in which a scalar Lagrange multiplier forces 
the d'Alembertian of a second scalar give R. That version of the model has 
additional degrees of freedom that the original proposal lacks [21], so it is 
important to examine the reconstruction problem in both formulations. 

This paper is organized as follows. In section [2] we summarize the model 
of Ref. [T2] and its specialization to cosmology. In section [3] we outline 
the steps of the reconstruction of the function f(X) once a given cosmology 
is chosen. In section HJ we show how to obtain the function / suitable to 
reproduce the ACDM cosmological evolution with the same matter content 
but no cosmological constant. We then summarize our results and conclude 
in section O 




(10) 
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2 Nonlocal Cosmology 



The modified gravity Lagrangian (fTUj) has already been given. It remains to 
state that □ is the covariant scalar d'Alembertian, 



By its inverse we mean the retarded Green's function. For the FLRW geom- 
etry (J4j) the action of this inverse on some function of time W(t) takes the 
simple form, 



1 

n 



w 



(t) 



o a A [u) Jo 



dva 3 (v)W(v) . 



(12) 



The metric g^ u is assumed to be minimally coupled to matter. 

The model is actually defined by its field equations, which are obtained by 
varying the gravity and matter actions with respect to the metric and then 
employing the partial integration trick explained in [TT]. This produces causal 
and conserved field equations, like those of the more rigorous Schwinger- 
Keldysh formalism [22] . These equations take the form, 



G>„ + AG 



N 1 fjLU-, 



(13) 



where is the matter total energy momentum tensor (including a possibly 
non vanishing cosmological constant), and AG Mi , comes from varying the non 
local / term in the above action (fTUl) . 

In the following, we will restrict the metric to be of FLRW form (jlj). With 
this ansatz, the field equations f|T3|) take the form 



3H 2 + AG, 



oo 



-2H - 3H 2 + — -5 ij AG. 
3a 2 



''.J 



8nGp, 
8ttGP, 



(14) 
(15) 



p and P being respectively the matter total energy density and pressure, and 
the non local pieces reading 



AG, 



oo 



3H 2 + 3Hd t \f(—R 



1 

□ 



1 



-R 



Rf 

+ k(i*)x*(i 



(16) 



4 




where /' denotes the derivative of / with respect to its argument. As already 
stressed, the left hand side of equation (fl3l) is conserved, and hence the 
first Friedmann equation (114p and the matter energy-momentum conservation 
equation, 

p + 3H(P + p) = 0, (18) 

are enough to ensure that equation ( fl5l) is fulfilled, as is the case with stan- 
dard Friedmann equations. 



3 General Reconstruction Technique 

We first note that the difference between field equations ( Tbfj) and ( fl5l) leads 
to a simple second order ODE for the function F defined as 




(19) 



where X is defined as 



X = -R. (20) 
□ 

This ODE reads 

F + 5HF + (6H 2 + 2H) (F + 1) = 8irG N (p - P) . (21) 

If one then assumes the matter content of the Universe (specified here-above 
by its total energy density p and pressure P) and its cosmological evolution 
(specified by some scale factor a(t)) to be chosen, the first step to recon- 
struct / is to solve Eq. (T2TT) . which allows to obtain F as a function of the 
cosmological time t. Then we invert relation (Tl9l) rewritten as 

■ dX 

UF = nf + Rf—, (22) 
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and yielding thus an ODE provided X is known as a function of t via equation 
( I20p . This allows to obtain / as a function^ of t. The last step is to invert 
X(t) to obtain t as a function of X. This allows to get f(X), the function 
which appears in the original action ffTUl) . 

In this process, some care has to be taken about the boundary conditions. 
First, we emphasize that our choice of the retarded Green function in the 
definition of the inverse of the d'Alembertian operator f[T2"l) does not permit 
inclusion of the extra zero modes that cause the ambiguities underlined in 
Ref . [21] about the local version of the model considered in [T8l [20] . More- 
over, one should make sure that whatever boundary conditions get imposed 
to integrate the ODEs are compatible with both equations <HM and (115]) . 
The need for this might seem confusing in view of the close relation between 
the two equations implied by conservation, 

— + 3H] (Eqn. [U) = -3H x (Eqn. Q22) . (23) 

However, it will be noted that relation (123!) involves a derivative of equation 
(JHJ). Had we based the reconstruction technique solely upon equation (fill) 
then equation (fT51) would have followed automatically. But our reconstruc- 
tion procedure instead employs the difference of equations (fT4"l) and (fT5]l . and 
this difference only defines equation (fllj) up to an integration constant. 

In the next section, we carry out these steps for the special case where the 
matter content of the Universe only consists of non relativistic and relativistic 
matter without any cosmo logical constant, while the cosmo logical evolution 
is that given by usual Friedmann equation with the same matter content plus 
a non vanishing cosmological constant, i.e. the one of the ACDM model. 

4 Specialization to ACDM 

We want to reproduce ACDM cosmology with the same matter content but a 
vanishing cosmological constant. Hence we assume that the Hubble param- 
eter H, appearing in equation (l2"Ti) . is a solution of the standard Friedmann 
equations with a cosmological constant and the same matter content as in 
our Universe: 

3H 2 -A = 8nGp, (24) 
2H + 3H 2 -A = -8nGP. (25) 

1 Denoting here by the same letter the function f(t) = f(X(t)) and the function f(X) 
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It is then easy to see that Eq (1211) simplifies dramatically to read now 



F + 5HF + (QH 2 + 2H)F = -6# 2 ft A , (26) 



where H is the Hubble parameter today, and Q A is defined as usual in terms 
of the cosmological constant A entering into equations (I24H25I) . i.e. as 

A 

m 



= ^H- (27) 



Notice that the matter energy density and pressure appearing on the right 
hand side of equation fl2T|) have cancelled against 6H 2 + 2H appearing on 
the left hand side. 

In the rest of the paper, we will further simplify the problem by consider- 
ing that the only matter content of the Universe is a sum of two components, 
one of non relativistic matter (with Q parameter Q m ) and one of relativistic 
matter (with Q parameter Q r ). It will then turn out to be convenient to use 
equations (I24j) and (|25|) to reexpress H and its time derivatives in terms of 
and Q m and to use the variable ( instead of t, defined in terms of the 
redshift z as 

C = l + z = -L (28) 
a{t) 

We also introduce the dimensionless Hubble parameter h(() given by 



h(C) = ^/Q A + Q m ( 3 + Q r (\ = H/H . (29) 
Equation (!26l) can be readily integrated^] to yield F in the form 

F = C 2 $ 2 , (30) 

where $ is given by 

and the integration constants $ e q and $g q are defined respectively as the 
values of $ and $' (here and in the following, a prime denotes a derivative 
with respect to () at the matter-radiation equality ( eq . In order to have a 



2 Using that the fact that C 2 is a homogeneous solution. 
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well behaved F at early times, we demand that F(Q goes to zero at large (. 



This fixes <3> eq and $' to be 



cq 



roo I roc I 

-60a/ ^Cl 777T / ^274^77- 



60/ 



cq 



00 1 



h e q J(cq ' ClMCO 



(31) 
(32) 



The resulting expression for $(C) is then given by 



roo I 

$(C) = -60 A / dCi 



00 I 
^2 



/c * KCi)^ C2MC2) 

and the large ( expansion of $(C) is, 



(33) 



*(C) 



50 r C 6 







(34) 



Any choice other than (13 lti32j) will result in a function F(Q actually growing 
for large (. 

Having obtained F(() we now turn to invert equation ffTOl) . This equation 
reads, 



C 2 /i 2 



d 2 



df 



d( 



dF 



(20 A + ±O m C 3 ) dC 



dC 



(35) 



Knowing F from (130]) it can easily be integrated once to yield the general 
solution 



-60, 





heqleq 


eq 


>>eq 






J 









C1 5 



• (36) 



This expression contains a new integration constant, (f') e q, and the function 
defined by 



- (120 A + 30 m Cf 
1 



(37) 
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Now, one also has from equation (1191) that 

F(n r(x(n) fV c? fV ( 12 "A+30nC2 3 ) x df 

" A^)) - / dCx wTT / ( ^ 6/( , x — . (38) 



where we have used the the definition of □ 1 given in equation (j!2p . In 
particular, we have that 



1 



R 



3^ 



C 2 (120 A +30 m C| 



-/ ^it^t/ ^ 2 V ; ,^r =X(Q. (39) 

■>c ft(Ci) ^ci C2MC2) 



which was used to obtain equation ( 1381) . The integral in Eq. ( |38l) will only 
make sense provided df/d( falls off faster that l/£ 2 for large (. This fixes 
the new integration constant (/')eq to be 



U'U = 2$ 



(:(/ 



^Ur dCl J£> _ 2 n Cl , (40) 



h eq Ieq \ J Ce 9 ClMCl) ^ fe 9 Cl 

and results in the new expression for df / d( reading 



00 



;i2Q A +3fi m Ci 3 )$(Ci 



2 . (41) 



C: 



If one integrates once more, one obtains another integration constant, namely 
the value of /(£) at C — Ceq- This constant can be fixed demanding that 
relation fl38|) hold at £ = ( eq . We get 



CO 



(I too (12fi A +3fi m C|) d/ 



uu - c*. - y^c, ^ y £i * a/(c 2 r x t • (42) 

Note that this relation is not self-referential for (f) eq because the factor of 
df /d(" on the right hand side does not involve (f) eq - With such a choice for 
(f)eq, the function f(x(() \ vanishes at early times, that is, for large (. We 
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can therefore construct / by simply integrating equation (j4ip to obtain 

df 



/(*<0) = -p&£, 

POO /°2 yCO 

+ 2 / ^1 777^ / ^C: 



C 2 4 MC 2 ) 



- (12fi A +3fi m C|)$(C 2 ) 

,2 



CI 



/•oo ^2 ^-2 

^2 



oo 1 



x AC 2 ) 

2 



C 2 4 MC 2 ) 



-36fi A / dCi T77T777T / 



CJ 



Ci 4 

o 



(43) 



Introducing the variable a = Ceq/C-, as we H as the parameter uo given by 
a; = Q^Qr/Q^, the above expression for / can be given in term of the two 
elliptic integrals J(a) and 1(a) defined by 



J(a) 
1(a) 



aai 

^1 + ai+uaf 



da i 



ai+4uaf 



1 + ai+uaf 3fi 



A3/2 

^ /(C). 



1/2- 



(44) 
(45) 



If we then define fi, f'2 and ^3 as the three contribution to / appearing in 
the right hand side of equation (H3|) . such that /(C) = /i(C) + /2(C) + /3(C)) 
one has 



/i(*(0) 



6a; / aai — ) =- 



(46) 



/■a 1 

-6a; / aai , 

-ofi+waf I (011) 



a.i\Jl + ai+ujai I(ai) Jo 

rot I 

12a; / aai , 

Jo a?A/l+a 1 4 



%, ^2) ,,(47) 



l+a2+a;a2 



(48) 
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X 



Q i (ai+co>a^ — a 2 — cua!) J (02) 
da>2 







(49) 



and one can use these expressions to numerically evolve / (X(Q). 

Having thus obtained / (X(Q) as a function of (, the last step is to get ( 
as function of X. Before doing so, it is also of interest to check that nothing 
goes wrong at late times. For that purpose we need the following small ( 



expansions, 



HO 

no 



4x701 3a m 



C 3 



A 



C 3 + 0(C 4 ), 

ln(C) + 0(l) 



-^5 + 0(l) 



(50) 

(51) 
(52) 



Now substitute these into each of the three terms in expression (HTj) for df/d( 
to get 



d f 2 



4Q 



8Q 



+ O(CH0) + (53) 







This implies that / approaches a constant at late times. In the above equa- 
tion (1431) . the three contributions appearing in the right hand side can be 
expressed using elliptic integrals. 

Let us now turn to obtaining an expression for X(Q. This also involves 
an elliptic integral. Indeed, equation (TSUI) reads 



X = - 



d( 2 



-i?(C 2 



C KCi) 



(54) 



where I is defined as in equation (|37|) . For a chosen set of parameters 
{Qa, Q m , fl r }, numerical evaluations of the right hand sides of equation (1541) 
and (H3|) can easily be obtained, from which one can get / as a function of X. 
The result is plotted in figure □ for {fi A ,fi m ,fi r } = {0.72,0.28,8.5 x 10~ 5 } 
which correspond to the latest WMAP values [2]. 

A simple analytic parameterization / an of the found function / is given 

by 



/anPO = 0.245 tanh(0.350F + 0.032F 2 + 0.003F 3 ) - 1 



(55) 
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f(X) 




Figure 1: Plot (solid blue curve) of the reconstructed function f(X) for the 
non local cosmology reproducing ACDM background cosmological evolution, 
with the same matter content but no cosmological constant. The parameters 
corresponding to the background cosmology are those of the latest WMAP 
release [2]. Circles indicate values of the function f(X) with the correspond- 
ing value of the redshift z indicated above. 
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where Y is defined by Y = (X + 16.5). It would hardly be distinguishable 
from the numerical solution should it be plotted together with the latter 
on figure [TJ These numerical and analytic expressions for f(X) allow us to 
resolve one of the major open questions about this class of models: do they 
make significant corrections to general relativity when expanded around flat 
space? The answer is "no." One can see from Figure 1 that the curve is 
almost flat near X = 0. From the analytic expression flBl)]) we compute, 



fL(o) 



.245 



0.350 + 0.064F + 0.009F C 



x sech z 



0.350F n + 0.032F 2 + 0.003F, 



o 



(56) 

10- 24 , (57) 
where Y = 16.5. So we find an utterly negligible linear correction. 



5 Discussion 



In this work we have presented a general method to reproduce a given arbi- 
trary cosmological evolution from a distorted non local form of the action for 
gravity as presented in Ref. [12]. This method was applied here to ACDM 
cosmology and we obtained the distortion function / that leads, via action 
( TTUj) . to exactly the same cosmological evolution as in ACDM, with the same 
matter content but no cosmological constant. It is very interesting to note 
that the function we obtain numerically is almost indistinguishable from a 
simple analytic form (|55|) . To be sure, this function contains some free pa- 
rameters — for example, the value X = —16.5 where the tanh passes through 
zero, or the scaling of its full variation by 0.49. However, these are all dimen- 
sionless numbers of order one. This is a consequence of two crucial properties 
of nonlocal models of type ([3]): 

• The onset of late time acceleration is triggered by the very recent 
cosmological transition from R « during radiation domination to 
R ~ 1/t 2 during matter domination; and 

• Even after this transition the nonlocal operator ^R ~ —\n(t/t eq ) 
grows very slowly. 

Several very interesting questions are left for future work. First, as far as 
cosmology is concerned, a natural question to address is if the model which 
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gives the same background evolution as ACDM can be distinguished from 
ACDM by considering observables that contain information in addition to 
the background evolution. This requires in particular working out the theory 
of cosmological perturbations for the non local model (see [19]). To do so, 
a good starting point is in fact this work, using for example, the analytic 
form (when necessary) of the reconstructed distortion function. Other issues 
concern the various tests of gravity one can consider, in particular those in 
the solar system or those involving binary pulsars. It would be extremely 
interesting to apply those tests to the framework of Ref. [12], and a first 
investigation along these lines has already been carried out in Ref. [T9] . 
Note in particular that the way cosmic acceleration is produced in the model 
considered here, is via a strengthening of the Newton constant encoded in 
the non local function /. However, the strengthening of the Newton constant 
we mentioned, strictly speaking only applies to cosmological distances, and 
things would be radically different, hence requiring a completely different 
analysis, inside gravitationally bound objects such as a galaxy or a cluster. 
This raises various questions about the effects of the non local operator CU 1 
inside matter. 
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